The Molecular Bases of the Interaction between a Saponin from the Roots of Gypsophila paniculata L. and Model Lipid Membranes

In view of the possible medical applications of saponins, the molecular structure of a GOTCAB saponin from the roots of Gypsophila paniculata L. was determined by NMR. The biological activity of saponins may depend on the interaction with cell membranes. To obtain more insight in the mechanism of membrane-related saponin function, an experimental and theoretical study was conducted. Ternary lipid systems composed of sphingomyelin, 1-palmitoyl-2-oleoyl-sn-glycero-3-phosphocholine, and cholesterol were used as models of mammalian cell membranes. The membrane–saponin interaction was studied experimentally by monitoring surface pressure in the monomolecular films formed at the air–aqueous subphase interface. The behavior of GOTCAB saponin in a water box and model monolayer systems was characterized by molecular dynamics simulations. The results obtained showed that, in the systems used, cholesterol had a decisive effect on the interaction between GOTCAB and phosphocholine or sphingomyelin as well as on its location within the lipid film.


Introduction
It has recently been shown that saponins may have therapeutic applications [1]. Cytostatic and cytotoxic effects of saponins in the case of malignant tumor cells have been described [2,3]. Other studies have shown that saponins exhibit adjuvant-active properties in vaccines as immunostimulatory complexes [4]. Based on these observations, an open cage-like immunostimulating complex of cholesterol, lipid, immunogen, and saponins from the Quillaja saponaria Molina bark (soap bark tree) was used as an active adjuvant in vaccines [5].
An example of a medical application of saponins is Saponinum album (SA; Merck), which is a crude mixture of triterpenoid saponins from Gypsophila paniculata L. and G. arrostii [6]. SA has been shown to enhance the cytotoxicity of type I ribosome-inactivating protein (RIP-I) saporin from Saponaria officinalis L. by 100,000-fold [7]. This interaction between saporin and SA is called "synergistic cytotoxicity principle". The combination of saporinbased targeted toxins and SA is used in anticancer therapy [8]. If a RIP-I is coupled with specific ligand that targets cancer-associated antigens, the cytotoxicity can even be increased by up to 4,000,000-fold in cell cultures [9]. This represents an important improvement in tumor therapy with targeted toxins as the side effects and costs of such a therapy could be lowered [8]. Böttger et al. [10] analyzed the structure-activity relationship of 56 saponins and their synergistic cytotoxicity with saporin and defined the concept of an ideal saponin. According to this concept, the ideal RIP-I synergistic saponin consists of an oleanane-type aglycone (gypsogenin or quillaic acid); a branched trisaccharide at C-3 consisting of β-D-glucuronic acid, β-D-galactopyranose, and β-D-xylopyranose; and a branched tetrasaccharide at C-28 consisting of deoxy sugars, such as β-D-fucose and/or α-L-rhamnose and acetyl residues.
The Gypsophila paniculata L. roots are rich in triterpenoid saponins. They contain about 4% dry weight in the fourth year of ontogenesis of the plant [11]. Phytochemical studies have shown that G. paniculata roots and SA comprise glucuronide oleanane-type triterpenoid carboxylic acid 3,28-bidesmosides (GOTCAB saponins) and monodesmosides [12][13][14][15][16][17]. The structure of GOTCAB saponins from G. paniculata and G. arrostii roots was first reported by Frechet et al. [12]. It was shown that gypsogenin and quillaic acid in the form of aglycones substituted at C-3 with a trisaccharide and with an oligosaccharide via fucose residue at C-28 are present in these molecules. It should be noted that the plant materials used by Frechet et al. were a mixture of both G. paniculata and G. arrostii. Unfortunately, SA cannot be used for further studies because Merck ceased its production in the 1990s. For the development of a combinatorial anticancer therapy in humans, saponins have to be isolated directly from the roots of G. paniculata. For the past 20 years, Gypsophila saponins have been purified by chromatographic techniques. A simple method for isolation of Gypsophila saponin-1641 from older batches of SA was proposed by Weng et al. [16]. Later, GOTCAB saponins were isolated using chemical degradation of the plant raw material followed by dialysis and HPLC purification [17]. Thakur et al. used an electrophoretic method to isolate saponins from SA [1]. Electrophoresis was also used to separate immune adjuvant saponins from Quillaja saponaria [18].
Different studies conducted by our group concerning Gypsophila species [19][20][21] have been previously published. Here, we present the results of a work aimed at elucidating the structure of a GOTCAB saponin from G. paniculata roots. This structural study was performed using NMR analysis. The molecular structure obtained was a starting point for further work concerning the role of saponins in pharmaceutical preparations. Because the interaction with cell membranes may play a role in the therapeutic activity of saponins, an experimental study was conducted using model mammalian cell membranes and molecular modeling.

NMR Analysis
1D and 2D NMR techniques were performed in order to elucidate the chemical structure of the isolated GOTCAB saponin (Supplementary Materials, Figure S2). The proton and carbon resonances of the aglycone moiety were assigned from the analysis of 1 H-1 H COSY, HSQC and HMBC spectra, the aglycone was identified as quillaic acid by comparison with literature data [14]. The 1 H NMR spectrum showed signals ascribable to six tertiary methyl groups at δ 0.72, 0.78, 0.83, 1.20, 1.50 ppm. Also evident were signals of H-23 at 9.55 ppm due to the aldehyde function at C-4 and of H-16 at δ 4.85 ppm due to the presence of a hydroxyl group at C-16. Complete assignment of the aglycone moiety is given in Table 1 According to spin-spin couplings and chemical shifts [12,14], one β-fucose (Fuc, Fuc-H1 at 5.59 ppm), two α-arabinose (Ara, Ara-H1 at 4.89 ppm; Ara', Ara'-H1 at 4.70 ppm), two β-xylose (Xyl, Xyl-H1 at 5.04 ppm; Xyl', Xyl'-H1 at 4.85 ppm), one α-rhamnose (Rha, Rha-H1 at 5.82 ppm), one β-galactose (Gal, Gal-H1 at 5.16 ppm), and one β-glucuronic acid (GlcA, GlcA-H1 at 4.66 ppm) were identified (assignments in Table 2). Glycosidation shifts were observed for Fuc-C2 (δ 73.4 ppm), Ara-C4 (δ 77.8 ppm), Xyl'-C3 (δ 84.6 ppm), Rha-C4 (δ 82.7 ppm), GlcA-C2 (δ 76.6 ppm), and GlcA-C3 (δ 84.2 ppm). The cross-peak in the HMBC experiment between GlcA-H1 (δ 4.66 ppm) and C-3 of the aglycone (δ 84.4 ppm) showed that the glucuronic acid was linked to the aglycone at C-3 position. Chemical shifts of Fuc-H1 (δ 5.59 ppm) and Fuc-C1 (δ 93.7 ppm) as well as the HMBC cross-peak between Fuc-H 1 (δ 5.59 ppm) and C-28 of the aglycone (δ 175.8 ppm) suggested that this sugar was involved in an ester linkage with the C-28 carboxylic group of the aglycone moiety. Key HMBC correlations were used in order to determine the position and the sequence of the sugar residues: between the proton signal at δ 5.16 ppm (Gal-H1) and the carbon resonance at δ 76.6 ppm (GlcA-C2), between the proton signal at δ 5.04 ppm (Xyl-H1) and the carbon resonance at δ 84.2 ppm (GlcA-C3), between the proton signal at δ 5.82 ppm (Rha-H1) and the carbon resonance at δ 73.4 ppm (Fuc-C2), between the proton signal at δ 4.85 ppm (Xyl'-H1) and the carbon resonance at δ 82.7 ppm (Rha-C4), between the proton signal at δ 4.89 ppm (Ara-H1) and the carbon resonance at δ 84.6 ppm (Xyl'-C3), and between the proton signal at δ 4.70 ppm (Ara'-H1) and the carbon resonance at δ 77.8 ppm (Ara-C4). The saponin was finally identified as corresponding to the structure presented in Figure 1.   It should be noted that extensive analysis of the NMR spectra described above allowed us to detect the presence of another saponin in the sample but at a very low concentration compared to other molecule described here (more than 5 times less according to the NMR signals integrals). Due to this low concentration, it was possible to distinguish the NMR signals of the two molecules, but it was not possible to completely assign the second one. Nevertheless, it can be assumed than this molecule was very similar to the GOTCAB saponin but had a gypsogenin moiety instead of quillaic acid and a different sequence of eight sugar moieties.
The structure obtained in this work ( Figure 1) was a starting point for further study on the GOTCAB-lipid membrane interaction.

Compression Isotherms
The interactions between GOTCAB saponin and the model cell membranes were investigated experimentally using Langmuir film technique [22][23][24]. The monomolecular films were formed with ternary SM/POPC/CHOL lipid mixtures containing biologically relevant 10 or 30% of cholesterol [25,26]; films containing 50 mol% of cholesterol were used as well. The SM/POPC equimolar ratio was kept constant. The final composition of ternary SM/POPC/CHOL monolayers was 45/45/10, 35/35/30, and 25/25/50. Some saponins show a high affinity for cholesterol, which is often linked to their hemolytic activity [27][28][29]. However, it was recently shown that membrane activity of the ginsenoside Rh2 saponin was boosted by the interaction with sphingomyelin but was worsened by cholesterol [30]. The membrane models used in this study were composed of a saturated (SM) and unsaturated (POPC) phospholipid as well as cholesterol; the latter is found in the raft structures of biological membranes. Sphingomyelin and cholesterol are involved in the formation of a raft-like liquid-ordered (L o ) phase, while POPC is involved in the nonraft liquid-disordered (L d ) phase [31][32][33].
The monomolecular films formed with ternary SM/POPC/CHOL mixtures on the pure water or aqueous GOTCAB subphases were characterized using compression isotherms. The interactions of GOTCAB with the model membranes were monitored by measuring the surface pressure (Π) as a function of molecular area (A). As a reference, a cholesterol-free equimolar SM/POPC monolayer was used. The isotherms obtained on the pure water subphase are presented in Figure 2. With increasing concentration of CHOL in the film, the isotherms shifted to lower molecular areas and the film rigidity increased, as indicated by the compression modulus ( Figure 2, the inset). This outcome agrees with many results published in the literature concerning phospholipid films containing cholesterol [34][35][36][37][38]. In the case of the subphases containing GOTCAB, the isotherms shifted to higher molecular areas compared to pure water. This effect indicates the interaction between GOT-CAB and the monolayers. As in the case of the pure water subphase, the SM/POPC/CHOL and SM/POPC films had different properties as well ( Figure 3). Indeed, the isotherms of the films containing cholesterol shifted to lower molecular areas compared to the binary films. In general, the films containing cholesterol were more rigid compared to the binary films ( Figure 3 insets and Figure S3 of Supplementary Materials). However, the ternary SM/POPC/CHOL 50 mol% films spread on GOTCAB were less stable compared to pure water, as indicated by a lower surface pressure at the collapse (Π coll ) of the monolayer (Figure 3). This effect was particularly well seen in the case of 80 mg·L −1 GOTCAB.
The Π coll values of SM/POPC/CHOL 50 mol% were around 35 and 29 mN·m −1 on water and 80 mg·L −1 GOTCAB, respectively. These observations indicate that cholesterol modifies the molecular interaction between GOTCAB and film-forming lipids. It may be supposed that GOTCAB forms different molecular associations with lipids depending on the presence of cholesterol in the film.

Molecular Dynamics Simulations
To obtain a better understanding of the experimental measurements, molecular dynamics simulations were used [39][40][41][42], which allowed the properties of the model membranes and molecular interaction within the system containing GOTCAB molecules to be studied at an atomic level.  As can be seen, GOTCAB HBs were formed preferentially with water. Approximately one GOTCAB molecule was involved in 14 HBs of water molecules. On the other hand, the number of HBs among GOTCAB molecules was less than one. Detailed analysis of the data obtained showed that intramolecular HBs were more common. The latter indicates that HBs are not responsible for the clustering the GOTCAB molecules.
The formation of the GOTCAB cluster is summarized in Table 3. The results correspond to a 12 ns production run. The numbers given are equal to the average number of clusters of a given size N (GOT) n observed during the last runs, where (GOTCAB) n is a cluster built of n monomers, and n = 1, 2, . . . denotes monomers, dimers, and so on. The numbers in each row are normalized to the total number of saponin molecules in the simulation box (27 in systems I-III): ∑ n n· N (GOT) n = 27. For example, the average number of monomers, dimers, etc. in system I at a given time was equal to 7.39, 4.69, etc., respectively. The same was true for small clusters, namely dimers and trimers. Larger clusters, with the exception of the tetramer for system III, were observed temporarily. For example, 0.58 for system II and N GOTCAB = 7 means a cluster was observed for almost 7 ns (0.58 × 12 ns ≈ 7 ns). This does not mean, however, that it is also the lifetime of such a cluster. A single cluster of this size was stable for 3 ns during the production run. Clusters of at least 9 monomers were not observed. The obtained data show that GOTCAB is present in water in the form of monomers and small aggregates. It has to be mentioned that the results obtained are qualitative and show a tendency. The HB pattern shown in Figure 4 indicates a hydrophobic origin of clustering. The structure of GOTCAB saponin shows many CH 3 groups attached to different residues (rhamnose, fucose, and aglycone). It can be expected that in clusters, such groups will be closer to each other. Therefore, a conformation with closely packed nearby CH 3 substituents will be preferred in the monomer. Such a spatial arrangement will minimize unfavorable interaction with the solvent. Figure 5 shows the radial pair distribution function for aglycone CH 3 as well as rhamnose and fucose. The first maximum at approximately 4 Å is a clear indicator of this preferred arrangement of the methyl groups. This situation is illustrated for the GOTCAB dimer in Figure 6. Radial pair distribution function for methyl carbon atoms of the aglycone moiety and methyl carbon atoms of rhamnose and fucose residues. The presented curve is an average over three independent simulations (system I-III, see Table 2).

Saponin Molecules at the Air-Water Interface
Partial density plots along the z axis, i.e., in the direction perpendicular to the monolayer surface, are shown in Figure 7. Panels (a) and (b) correspond to binary SM/POPC and ternary SM/POPC/CHOL 30 mol% monolayers in the presence of GOTCAB molecules, respectively.  Figure 7 illustrates densities of the terminal hydrophobic CH 3 groups in phospholipid lateral chains, hydrophilic PO 4 , N(CH 3 ) 3 , and OH functional groups. Moreover, the location of C=C double bonds is shown. Independently of the system, the outspread of the hydrophobic parts of the lipids towards vacuum was almost the same. The maximum of the partial density plot of terminal CH 3 groups of POPC (magenta area) and SM (violet line) in panel (a) were located in the same region, including the location of maxima. The same was true for the ternary system. An additional orange line representing the terminal CH 3 groups of CHOL was exactly in the same place as the remaining two lines. The differences were more pronounced in the hydrophilic region. Here, the partial density plots of PO 4 in POPC (orange area in both panels) shifted towards the water box interior compared to the SM PO 4 groups (red line). The hydrophobic part of CHOL, represented by the black line in panel (b), was located at the interface. The same behavior was found in POPC N(CH 3 ) 3 and SM N(CH 3 ) 3 groups, with the former pushed further towards the water bulk. A similar trend was observed in the SM/POPC and SM/POPC/CHOL monolayers formed on pure water (Supplementary Materials, Figure S4).
In Figure 7, the interaction of GOTCAB and SM/POPC or SM/POPC/CHOL is indicated by the relative positions of the corresponding maxima. It can be observed that the yellow areas (GOTCAB) presented in panels (a) and (b) are dissimilar. In the binary system (panel a), GOTCABs interacted with the hydrophilic and partially hydrophobic regions of the monolayer formed by different lipid moieties. In the ternary monolayer (b), the yellow area shifted towards the aqueous subphase, indicating that polar interactions between GOTCABs and film-forming lipids were favored. These results show that in the presence of cholesterol, GOTCABs do not penetrate into the monolayer.
The side view of the monolayers is presented in Figure 8. In panels (a) and (c), all molecules forming the monolayer are shown; in panels (b) and (d), POPC molecules are not displayed to facilitate visualization of the interaction between GOTCAB and the lipids. It can be observed that in the binary SM/POPC system, GOTCABs formed aggregates in the aqueous phase (see Figure 5), while SM-GOTCAB complexes were formed via polar and apolar interactions in the film (Figure 8a,b). In the case of GOTCABs immersed in water, the small number of hydrogen bonds was due to interactions between sugar OH groups and oxygen atoms in the phospholipid PO 4 groups. However, most HBs were formed with water molecules. In the SM-GOTCAB complexes, the aglycone moieties interacted with the SM hydrocarbon chains, while the sugar residues were present in the hydrophilic region of the SM/POPC film. The results obtained showed that -NH-and OH groups in SM were involved in the reorientation of GOTCABs in the monolayer ( Figure S5). Indeed, sugar and aglycone hydroxyl oxygens formed hydrogen bonds with nitrogen and oxygen atoms of the -NH-CO-group, respectively, stabilizing GOTCAB in the apolar region of the lipid film (Figures S5 and S6). The -O-CO-groups in POPC were practically not involved in hydrogen bond formation (compare both curves in Figure S5). It can be noted that the effects observed resemble the flip/flop mechanism described in the case of the CX1 calixarene derivative [40,43]. bonds with nitrogen and oxygen atoms of the -NH-CO-group, respectively, stabilizing GOTCAB in the apolar region of the lipid film ( Figures S5 and S6). The -O-CO-groups in POPC were practically not involved in hydrogen bond formation (compare both curves in Figure S5). It can be noted that the effects observed resemble the flip/flop mechanism described in the case of the CX1 calixarene derivative [40,43]. In the ternary SM/POPC/CHOL system, GOTCAB cluster distribution resembled that observed in pure water (compare Figure 5 and red curve of Figure S7). Here, GOTCABs interacted with the POPC polar heads and stayed at the lipid-aqueous interface (Figure 8c,d). The latter observation was in accordance with that based on the density profile analysis (Figure 7). Overall, the results obtained with modeling indicate that cholesterol hinders penetration of GOTCABs into the monolayer and makes formation of GOTCAB-SM complexes less likely. This observation is consistent with some results published in the literature showing that, under the conditions used, cholesterol is preferentially associated with sphingomyelin [44].

Materials
The saponin used in this study was purified as described previously [21]. The molecular structure of the GOTCAB saponin is shown in Figure 1. Sphingomyelin (Egg SM, ~99.9% pure) and 1-palmitoyl-2-oleoyl-sn-glycero-3-phosphocholine (POPC, ~99.9% purity) were from Avanti Polar Lipids. Cholesterol (CHOL) and chloroform (~99.9% purity) used for preparing phospholipid solutions were from Sigma-Aldrich. Pure water and aqueous GOTCAB solutions with concentrations of 8 and 80 mg L −1 , respectively, were used as subphases for the monolayer experiments. The subphases were prepared with MilliQ water that had a resistivity of 18.2 MΩ·cm at 25 °C, surface tension of 72.8 mN m −1 at 20 °C, and pH of 5.6. In the ternary SM/POPC/CHOL system, GOTCAB cluster distribution resembled that observed in pure water (compare Figure 5 and red curve of Figure S7). Here, GOT-CABs interacted with the POPC polar heads and stayed at the lipid-aqueous interface (Figure 8c,d). The latter observation was in accordance with that based on the density profile analysis (Figure 7). Overall, the results obtained with modeling indicate that cholesterol hinders penetration of GOTCABs into the monolayer and makes formation of GOTCAB-SM complexes less likely. This observation is consistent with some results published in the literature showing that, under the conditions used, cholesterol is preferentially associated with sphingomyelin [44].

Materials
The saponin used in this study was purified as described previously [21]. The molecular structure of the GOTCAB saponin is shown in Figure 1. Sphingomyelin (Egg SM, 99.9% pure) and 1-palmitoyl-2-oleoyl-sn-glycero-3-phosphocholine (POPC,~99.9% purity) were from Avanti Polar Lipids. Cholesterol (CHOL) and chloroform (~99.9% purity) used for preparing phospholipid solutions were from Sigma-Aldrich. Pure water and aqueous GOTCAB solutions with concentrations of 8 and 80 mg·L −1 , respectively, were used as subphases for the monolayer experiments. The subphases were prepared with MilliQ water that had a resistivity of 18.2 MΩ·cm at 25 • C, surface tension of 72.8 mN·m −1 at 20 • C, and pH of 5.6.

NMR Experiments
NMR experiments were performed on a Bruker Avance III spectrometer operating at 9.4 Tesla (400 Mhz and 100.6 MHz for 1 H and 13 C, respectively) using a Bruker 5 mm BBFO probe. Pulse widths were 14.1 and 10.5 µs for 1 H and 13 C, respectively. Samples were dissolved in 400 µL of 50% pyridine-d 5 and 50% D 2 O (pyridine-d 5 was used as reference, with highest field signals at 7.19 ppm for 1 H and 123.5 ppm for 13 C), and all experiments were performed at 313 K. Standard experiments were performed using the Bruker software package: 1 H, 13 C decoupled from proton, 13

Surface Pressure-Area Isotherms
The surface pressure (Π) measurement was carried out using a KSV 2000 Langmuir balance (KSV Instruments, Helsinki, Finland). A Teflon trough (58 × 15 × 1 cm) with two hydrophilic Delrin barriers providing a symmetric compression was used in all experiments. Surface pressure measurements were carried out using the Wilhelmy plate method. The apparatus was closed in a Plexiglas box, and the temperature was kept constant at 20 • C. All impurities were removed from the subphase surface by sweeping and suction. Monolayers were spread from chloroform solutions of accurate lipid concentrations using a microsyringe (Hamilton Co., Reno, NV, USA). SM/POPC mole fraction was 0.5 (50 mol%). CHOL mole fraction in the SM/POPC/CHOL mixtures was 0.1, 0.3, or 0.5 (10, 30, or 50 mol%, respectively).
After the equilibration time of 20 min, the films were compressed at the rate of 2.5 mm·min −1 by two symmetrically moving barriers. A PC computer and KSV software were used to control the experiments. Each compression isotherm was performed at least three times. The accuracy of the results was ±0.1 Å 2 for mean molecular area and ±0.01 mN·m −1 for surface pressure measurements.

Molecular Dynamics Simulations
The behavior of GOTCAB saponin in a water box and model monolayer systems was characterized by molecular dynamics (MD) simulations. The former was composed of 27 GOTCAB molecules distributed randomly on 3 × 3 × 3 grid. In the symmetric monolayer system, 9 GOTCAB molecules were placed 10 Å under/above the equilibrated monolayers on 3 × 3 grid. Two model monolayers were built: a binary monolayer comprising SM and POPC (SM/POPC monolayer) and a ternary monolayer with addition of CHOL molecules (SM/POPC/CHOL monolayer). The initial configurations for both types of simulations are shown in Figure 9. The GOTCAB saponin in the water box is shown in panels a and c. The ternary monolayer model with GOTCAB is shown in panels b and d. All-atom CHARMM force field was applied [48,49]. The missing parameters of GOTCAB saponin were derived using the force field toolkit of VMD graphical software [50,51]. The initial structure of GOTCAB, namely one of the many conformers, was optimized at the B3LYP/6-31G(d) level of theory. Due to its size, the saponin was divided into smaller fragments, each retaining the chemical identity of the immediate environment. A force field toolkit was applied to each fragment separately. Such a procedure is in line with the standard assumption of molecular dynamics sumulations, i.e., the transferability of force field parameters. The parameters derived in this way were combined together, including the charge distribution. The MD calculations were carried out using the NAMD package [52] with periodic boundary conditions imposed on the system. Langevin thermostat and barostat were used to control the temperature and pressure, respectively. van der Waals interactions were switched off at 12 Å. Electrostatic energy was calculated with the PME method [53]. A time step of 1 fs was used in all simulations. TIP3P model was adopted for water [54]. The VMD software was also applied to visualize the trajectories and to compute selected properties, including hydrogen bond characteristics. It was assumed that the distance from the hydrogen donor (HD) atom to the hydrogen acceptor (HA) atom was less than or equal to 3 Å and that the HD-H-HA angle was not greater than 20 degrees.
In all GOTCAB/H 2 O calculations, the number of atoms (N) in the simulation box was greater than 160,000. Three independent simulations, each of 120 ns, were performed for GOTCAB saponins. The systems were equilibrated in (N, p, T) ensembles for 60 ns with T = 293 and p = 1 bar. The second half of calculations were performed in canonical ensembles. The last 12 ns runs were used to compute the average quantities. The binary and ternary monolayers consisted of 100 lipid molecules. The number of atoms in all systems were greater than 120,000.
Monolayers were placed on two sides of a water slab (xy cut). After the initial 10 ns canonical ensemble simulations, the calculations were carried out at constant temperature (T = 293 K), normal pressure (p n = 1 bar), and surface tension. Two values of surface tension (γ = 50 and 30 mN·m −1 ) were considered. The (N, T, p n , γ) simulations were performed for 120 ns. These calculations were followed by 60 ns canonical ensemble calculations with the box size fixed from (N, T, p n , γ) calculations (average from the last 12 ns). Thus, the total simulation time for each system was 190 ns. The data presented in the paper were computed using the last 24 ns of the final canonical ensemble run being considered as the production run. The picture presenting the equilibration process is shown in Figure S1 in the Supplementary Materials (after 100 ns, the total potential energy fluctuated around the average value).

Conclusions
The results obtained in this work clarify the mechanism of the interaction between a pure GOTCAB saponin with defined molecular structure and model cell membranes containing phospholipids or phospholipids and cholesterol. In the monolayers prepared with the two phospholipids, namely sphingomyelin and phosphatidylcholine, GOTCAB established polar interaction with the phospholipid polar heads (PO 4 groups) via sugar moieties (OH groups), and apolar interaction between the aglycone moiety and hydrocarbon chains present in sphingomyelin. Obviously, the interaction with the latter was favored compared to the phosphatidylcholine chains. The conjunction of polar and apolar interactions allowed stabilization of the GOTCAB in the apolar part of the monolayer. In the ternary monolayers, the interaction between cholesterol and sphingomyelin was favored compared to GOTCAB-sphingomyelin. This effect may be explained by a better cholesterol match between the hydrocarbon chains compared to the saponin aglycone moiety. Consequently, in ternary films, GOTCAB was stabilized in the region of the lipid polar heads. It can be concluded that, in view of its effectiveness, competition with cell membrane cholesterol has to be considered when medical, saponin-based preparations are developed. The structure and purity of the saponin used may be crucial for the final biological effect of such preparations.

Funding:
The calculations were carried out at the Faculty of Chemistry of the Jagiellonian University and on the PL-Grid infrastructure at ACK CYFRONET with the support of the project "HPC Infrastructure for Grand Challenges of Science and Engineering". K.K. acknowledges the support of InterDokMed project no. POWR.03.02.00-00-I013/16.

Institutional Review Board Statement: Not applicable.
Informed Consent Statement: Not applicable.